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摘要 : 本 文 提 出 一 种 考虑 流体 热 物 性 变化 的 格子 - 玻 尔 效 曼 通 量 求解 器 。 该 求解 器 能 够 捕捉 流体 物性 参数 的 变化 对 流动 和 传 热 
的 影响 ， 在 保留 了 传统 格子 - 玻 尔 兹 曼 方法 固有 优势 的 同时 ， 克 服 了 在 网 格 划分 、 边 界 条 件 处 理 等 方面 的 缺陷 。 本 文采 用 该 求 
解 器 模拟 了 同心 环 型 腔 内 由 温差 驱动 的 自然 对 流 ， 讨 论 了 不 同 温差 条 件 下 的 流动 和 传 热 特性 ， 分 析 了 物性 变化 的 影响 。 结 果 
表明 : 常 物性 解 低估 了 换 热 设备 的 传 热 性 能 ， 而 且 Ra 数 越 大 ， 偏 差 越 显著 。 
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A Variable Property-based Lattice Boltzmann Flux Solver 


CAO Yu-Hui! 
(1. College of Sciences and Engineering, University of Chinese Academy of Sciences, Beijing 100049, China) 


Abstract: A variable property-based lattice Boltzmann flux solver is proposed in the present paper. This solver is 
capable of capturing the effects of the variation in fluid properties on flow and heat transfer characteristics. It retains 
inherent advantages of conventional lattice Boltzmann method and overcomes many drawbacks such as the difficulty 
in using non-uniform meshes, the inconvenience of dealing with boundary conditions, etc. In the present paper, the 
natural convection induced by a radial temperature difference in a horizontal concentric annulus is simulated by the 
solver, the flow and heat transfer characteristics obtained under different temperature difference conditions are 
studied, and variable property effects are discussed. It is found that the commonly-concerned constant property 
solution (CPS) underestimates the heat transfer performance of heat exchangers, and the deviation of the CPS from 
the variable property solution becomes more and more notable with the increase of the Rayleigh number Ra. 
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0 引言 展 性 ， 在 模拟 大 规模 复杂 非 线性 流动 问题 时 具有 独 
= 特 优势 ， 已 被 成 功用 于 处 理 各 类 复杂 流动 问题 上 
格子 - 玻 尔 兹 曼 方 法 (LBM) 基于 介 观 动 理 7 ， 如 工程 传 热 传 质 、 溃 流 、 多 组 分 、 多 相 流 、 界 

论 ， 其 中 的 演化 过 程 物理 图 像 清晰 ， 计 算 简单 ， 容 。 ” 面 动力 

易 编程 ， 计 算是 局 部 的 ， 具 有 良好 的 并 行 性 和 可 扩 
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学 、 化 学 反应 与 燃烧 等 。 

最 近 两 年 ，Shu 等 人 8 引 提 出 一 种 格子 - 玻 尔 兹 点 需要 强调 指出 。 第 一 ， 考 虑 了 流体 的 昔 性 、 导 热 
曼 通 量 求解 器 〈LBFS) ,并 将 其 推广 应 用 于 单 相 流 系数 、 比 热 容 以 及 体积 力 中 流体 密度 随 温度 的 变 
体 、 多 相 流体 的 流动 与 传 热 的 数值 研究 UC. 在 化 。 第 二 ， 可 以 很 方便 地 采用 非 均匀 网 格 。 第 三 ， 
LBFS 中 ， 通 过 对 格子 - 玻 尔 效 曼 方程 进行 C-E 分 析 存储 量 小 ， 只 需 存 储 平 衡 分 布 函数 。 第 四 ， 普 适 性 
得 到 守恒 的 微分 方程 ， 然 后 采用 有 限 体积 方法 对 微 好 ， 方 程 中 的 体积 力 以 及 由 物性 变化 引入 的 附加 项 
分 方程 进行 离散 化 处 理 ， 单 元 界面 上 的 数值 通 量 利 不 会 影响 数值 通 量 的 估算 过 程 。 近 期 ， 本 文 作者 提 


用 格子 - 玻 尔 效 曼 方程 的 局 部 解 进行 估算 。 有 关 出 一 种 可 以 同时 考虑 流体 密度 变化 的 算法 
LBFS 的 详细 介绍 ， 可 以 参考 文献 [81。 (VPLBFS)， 在 文献 1 中 做 了 详细 介绍 和 系统 检 


变 物性 的 影响 是 过 去 几 十 年 理论 研究 和 数值 研 验 。 


究 的 热点 问题 之 一 。 已 有 研究 结果 表明 : 流体 热 物 ed E SUE 
性 随 温度 的 变化 对 流动 不 稳定 性 有 比较 显著 的 影响 。 工 变 物性 格子 - 玻 尔 兹 曼 通 量 求解 器 


H3-16]， 同 时 对 传 热 性 能 的 影响 也 不 容 忽视 0 241 。 1.1 宏观 控制 方程 
然而 ， 采 用 LBM 研究 这 类 问题 的 工作 525-26 比 较 少 对 于 变 物性 流体 的 流动 和 传 热 问 题 ， 在 忽略 理性 热 
见 ， 而 且 都 只 能 采用 均匀 网 格 。 耗 散 和 压缩 功 的 条 件 下 ， 控 制 方程 可 以 写成 如 下 形 
本 文 将 LBFS 推广 用 于 变 物性 流体 流动 与 传 热 式 : 
过 程 的 数值 研究 。 关 于 本 文中 的 变 物性 LBFS， 有 四 
Op+V:(pu)=0 (1) 
ó (pu) * V -(puu) = -Vp  V| (Vu «(Vuy)|«F (2) 
1 
— —(8,T - u- VT) - V-(xVT) (3) 
zi ) 


1.2 标准 LBE 和 多 尺度 C-E 分 析 
类 似 于 常 物性 流体 ， 变 物性 流体 流动 的 标准 LBE 也 
和 比热容 C, 都 是 温度 的 函数 ， 是 体积 力 ， 只 有 可 以 写成 方程 (4) 的 形式 : 


体积 力 中 出 现 的 流体 密度 是 温度 的 函数 。 


在 本 文 的 研究 工作 中 ， 流 体 的 黏度 UL ， 导 热 系数 


Jfa (r, DL. t) 


f 


a=0,1,......,N (4) 


falr +e, E+ O) = fats t) + 


其 中 ，7 是 坐标 位 置 ， 0 是 粒子 运动 的 时 间 步 长 ， 


2 2 
£56.) pa, | 14 26, Uu) "GI | quy 
C; 2c, 


ROUIET FIRES, e, 是 格子 模型 中 心 方 


向 上 的 粒子 速度 ， X 是 粒子 速度 的 个 数 ，z， 是 弛 tt, o, 是 加 权 系数 ，c EER, o, c, 的 数值 


取决 于 所 选用 的 格子 模型 。 对 于 二 维 问题 ， 最 常用 

的 格子 模型 是 D2q9 模型 ， 图 1 给 出 了 D2q9 模型 的 

应 的 平衡 分 布 函 数 ， 具 有 方程 (5〉 的 形式 : 示意 图 ， 各 个 方向 上 粒子 的 速度 具有 C6 的 形 
式 : 


Hel], fa 方向 上 的 密度 分 布 函数 ， je 是 对 
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ee 二 


a 


(cos[(a —1)z / 2], sin[(e — Dz / 2] )c, 
V2 (cos[(a —5)z / 2 z / 4], sin[(a —5)z / 2 z / 4] c, 


c, - cl NB ， 加 权 系数 w Ha, = AIO , 
Q,-0,-0,-0,-1/9, 


f, TH fA ESTEE (7) 和 CD : 


对 方程 (4) fH Taylor 展开 ， 得 到 下 面 的 方程 (9) 


a 


图 1 D2q9 格子 模型 示意 图 


Fig. 1 Schematic of the D2Q9 lattice model 


在 多 尺度 C-E 分 析 中 ， 分 布 函数 大 以 及 大 对 
时 间 和 空间 的 偏 导 数 按照 以 下 形式 进行 多 尺度 展 


开 ， 

f," ref +e fo (10a) 
Sh are A (10b) 
a at, t 

V= a (10c) 
其 中 ，& 是 一 个 小 量 ， 其 大 小 正比 于 努 森 数 。 

将 方程 (10) 代入 方程 (9)， 经 过 整理 得 到 下 面 
的 方程 : 


2 
Fax -V 
2 \ ot 
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a=0 

a =1,2,3,4 

a=5,6,7,8 
O; = 0, = 0, = 04 =1/36 。 
p= Dif 232 


1 dl 2| 
forget )eo(ór)-o 


op 
E MS =0 
ap +V- (eu) 


27 + Y-(puu)=-Vp+V/ wus 


其 中 ， 


jsp (z, -0.5)6 


对 于 温度 场 ， 引 入 分 布 函数 g ， 


a 


(6) 


(7) 


(8) 


(9) 


(11) 


(12) 


(vu) | 


(13) 


， 它 的 平衡 分 布 函 


数 记 作 g“。 在 本 文中 ， 对 温度 分 布 函 数 g, 采 用 了 


D2q4 模型 ， 对 g, 满足 的 LBE 


C-E 分 析 之 后 ， 可 


方程 (15), 


<+V-(uP)=V-(DVT) 


1 
D- (n, -0.55 


1.3 变 物 性 LBFS 的 控制 方程 
从 标准 LBE 出 发 得 到 的 宏观 控 


以 得 到 等 恒 律 


进行 Taylor 展开 和 


(14)， 宏 观 控制 


TIVI TER PRESE TR] r。 之 间 的 关系 式 


(14) 


(15) 


(16) 


HAE C110. C12) 


Al (150 跟 变 物性 流体 自然 对 流 的 控制 方程 《1)、 


(2) 和 (30 是 有 差别 的 。 为 了 恢复 方程 (DD - 
3)， 同 时 又 能 保留 格子 - 玻 尔 兹 曼 模 型 和 守恒 律 
的 简洁 形式 ， 本 文 将 由 标准 LBE 得 到 的 守恒 律 和 关 
系 式 直接 代入 宏观 控制 方程 COD - (3) ,从 而 得 到 
下 面 的 方程 : 


(17c) 
对 于 二 维 问题 ， 在 变 物性 LBFS 中 ,方程 (17) 
可 以 改写 成 以 下 形式 : 


OW 0E OF je 
Ot Or oy 

(18) 

其 中 ， 


W={p, pu, pv, TY ; 

(19a) 

Esp m dp 

C19b) 

F={P, T, n, 9]; (19¢) 
D T 

Fé ={o 0, F, -wrv2| 
K 

(19d) 


N 
=, eq 
P. = X End, a 
a=0 


(20a) 
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a=0 L 277 
II ze e Fr 1- dud 
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(20d) 
ONE k s (20e) 
II = Jg bae ipse] goa 
y Deasa Se zc | 
1 (20f) 


在 变 物性 LBFS 中 ， 采 用 有 限 体积 方法 求解 方程 
(18)， 其 中 ， 每 个 单元 中 心 的 宏观 物理 量 通过 直 
接 求解 离散 化 方程 得 到 ， 相 邻 单元 界面 上 的 数值 通 
量 利用 LBM 的 解 通过 局 部 重 构 进行 估算 。 按 照 有 限 
体积 方法 ， 对 方程 (18) 在 控制 体 Q2, 上 进行 积分 ， 
得 到 变 物性 LBFS 需要 求解 的 控制 方程 (21)。 
dW, 1 
(2D 

其 中 ，AV 是 控制 体 Q, 的 体积 ，AS, 是 第 大 个 交界 
H R) 的 面积 (长度)。R 是 第 上 个 交界 面 
CX) 上 的 数值 通 量 。 估 算 Ri 的 具体 做 法 ， 请 参考 


4 阶 龙 格 - 库 塔 方法 得 


SHR 8]. AIR, 之 后 ， 采 
到 方程 (21) 的 数值 解 。 
2 变 物 性 LBFS 算法 检验 


为 了 检验 变 物性 LBFS 的 可 靠 性 ， 在 模拟 同心 
环 型 腔 内 的 自然 对 流 之 前 ， 首 先 重 复 了 文献 Cun 
的 工作 ， 模 拟 了 方 腔 内 的 自然 对 流 。 空 腔 的 上 、 下 
壁面 是 绝热 条 件 ， 左 、 石 两 侧 分 别 是 温度 恒定 的 高 


温 壁 面 和 低温 壁面 ， 温 度 分 别 记 作 开 MT. H 


` 


(22) 定义 了 一 个 无 量 纲 温 差 比 2 ， 用 于 衡量 温差 
的 大 小 。 


0 - (T, - T, )/T, (22) 


在 本 文 的 三 个 算 例 中 ， 瑞 利 数 Ra tH 10°, 37 
朗 特 数 Pr EX 0.71, T, EX 300K, O 的 取 值 分 别 
为 0.0101, 0.403 和 0.8, FFL, T, 分 别 为 


303.03K, 420.9K 和 5$40K。 流 体 的 热 物 性 参数 ， 黏 


度 4L， 比 热 容 C, 和 导热 系数 k ， 都 是 温度 的 函 


数 ， 函 数 表达 式 采用 文献 
流体 的 密度 仍 采 月 


中 的 形式 ， 体 积 力 中 


H Boussinesq 假设 。 


图 2， 不 同 温差 条 件 下 稳 态 解 的 比较 。 瑞 利 数 Ra =10° 。 
速度 无 滑 移 条 件 ， 左 、 右 侧 的 恒温 条 件 和 上 、 下 壁 
面 的 绝热 条 件 。 初 始 时 刻 ， 流 体 处 于 静止 状态 ， 流 
体 的 温度 等 于 低温 壁面 温度 。 
展示 了 不 同 温差 条 件 下 稳 态 解 之 间 的 比 

较 ， 包 括 流 线 和 等 温 线 。 如 图 2 所 示 ， 流 体 热 物性 


图 2 


的 变化 对 流 场 和 温度 场 都 有 比较 明 
温差 的 增 大 ， 变 物性 
分 析 流 体 热 物性 的 变化 对 传 热 的 


高 温 壁面 


的 影响 逐渐 增强 。 


显 b 影响 ， 随 着 
为 了 定量 地 


区 啊 ， 计 算 了 左 侧 


上 的 平均 努 塞 尔 数 Nu 。 


图 3 给 出 了 不 同 


温差 条 件 下 得 到 的 Nu 随时 间 的 变化 。 如 图 3 所 
示 ， 温 差 为 0.0101 时 的 Nu 非常 接近 常 物性 解 中 的 


数值 ， 随 着 温差 的 增 大 ， 传 热 速率 逐渐 增 大 。 图 2 


和 图 3 中 展示 的 结果 与 文献 2 中 的 结果 吻合 ， 
了 效 捕 捉 到 流体 热 物性 随 温度 的 变 


变 物 性 LBFS 能 


化 对 流 场 和 温度 场 造成 的 影响 。 
3 同心 环 型 腔 内 的 自然 对 流 


表明 


本 文采 用 变 物 性 LBFS 模拟 了 同心 环 型 腔 内 由 径 向 
温差 驱动 的 流动 和 传 热 过 程 。 同 心 环 的 内 环 半径 和 


外 环 半 径 分 别 记 作 R, 和 R,， 纵 横 比 

rr=R,/ R 22.6. Ay P RETI Bice T, M T, o 
跟 第 2 部 分 相同 ， 通 过 不 同 温差 条 件 下 数值 结果 之 
间 的 比较 讨论 流体 热 物性 变化 的 影响 。 工 恒 为 
300K， 无 量 纲 温 差 比 9 分 别 为 0.0101，0.403 和 
0.8。 普 朗 特 数 Pr 恒 为 0.71。 初 始 时 刻 流体 处 于 静 
止 状态 ， 流 体温 度 为 也 。 速 度 边 界 采 用 无 滑 移 条 


fr. i8 


径 向 采 上 月 


非 均匀 网 格 ， 沿 周 向 采 


ISHAM. 


EW: 流 线 ， 右 侧 : 等 温 线 。 
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— 6,=0.0101; 


-- 6,-040. -.- 6,-08. 


Fig. 2 Comparison of steady-state solutions obtained under 


different temperature difference conditions. Ra = 10^. 
Left: streamline, Right: isotherm. —— 6, 20.0101; 


0,-0405; =- 86,-08. 
跟 传 统 的 LBM 不 同 ， 变 物性 LBFS 可 以 很 方便 
地 采用 非 均 匀 网 格 进行 数值 求解 ， 而 且 不 必 通 过 分 
布 函数 给 定 边 界 条 件 ， 可 以 直接 利用 宏观 物理 量 给 
出 。 对 于 上 述 条 件 下 空 腔 内 的 自然 对 流 ， 采 用 变 物 


性 LBFS 进行 求解 ， 经 网 格 独立 性 检验 ， 发 现 使 用 
50x 50 的 非 均匀 网 格 足以 得 出 准确 的 稳 态 解 。 另 
外 ， 对 于 边界 条 件 的 处 理 ， 直 接 给 定 四 个 壁面 上 的 


3.1 常 物性 解 


首先 ， 采 用 变 物 性 LBFS 的 简化 


TN, MRA 


度 、 导 热 系 数 、 比 热 容 随 温度 的 变化 ， 得 到 了 同心 


7 


环 型 腔 内 自然 对 流 的 常 物性 解 。 网 格 数 为 180x 40 


。 图 4 给 出 了 常 物性 解 中 速度 分 量 和 等 温 线 的 分 
布 。 瑞 利 数 Ra = 50000 ， 流 体 的 物性 参数 ， 如 : 


黏度 4L， 比 热 容 C, 和 导热 系数 上 ， 取 温度 为 ZT. 时 


的 值 。 这 些 定性 的 结果 跟 文献 9 中 一 致 。 在 本 文 
中 ， 把 常 物性 解 作为 参考 ， 在 下 文中 将 通过 定量 的 
比较 ， 分 析 变 物性 的 影响 。 


— — 0 7 0.0101 (Present) 
9 0-0.0101 (Leal et al., 2000) 

----070.403 (Present) 
0—0.403 (Leal et al., 2000) 

(Present) 
(Leal et al., 2000) 


eM 0 - 0.8 
^ 0-08 


0.01 0.1 


Time (s) 


i 上 了 


图 3 左 侧 高 温 壁 


F 均 努 塞 尔 数 的 时 间 演 化 。 


Fig. 3 The temporal evolution of the average Nusselt number Nu 


calculated on the left hot wall. 
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4 常 物性 解 中 的 速度 分 量 和 等 温 线 。Ra = 50000 


Fig. 4 Distributions of velocity components and isotherms in constant property solution. Ra = 50000 


3. 2 变 物性 的 影响 响 逐 渐 增 强 。 
在 这 部 分 ， 采 用 变 物性 LBFS 模拟 了 三 个 不 同 Q 条 件 为 了 定量 地 分 析 物 性 变化 对 传 热 速率 的 影响 ， 

司 心 S ps E 然 对 ; PES, i , 热 系 Y [NO m ` i 
PRELDGMBISPRDHARADE. MENRE SRR TIR EER RS A, RAR 
数 和 比热容 都 是 温度 的 函数 ， 函 数 关系 仍 采用 文献 
[18] 中 的 形式 。 Ra 和 Pr BAA janes 和 0. Js [28], Kup 可 以 被 定义 为 (23) 的 形式 

5 展示 了 不 同 Q 条 件 下 数值 解 之 间 的 比较 ， 包 
HERRERA RAIN AT. BARRAR reer) pear TN 
似 ， 热 物性 的 变化 对 速度 场 和 温度 场 都 有 比较 明显 ” 2r0r-D gr 
影响。 而 且 ， 随 着 温差 比 的 增 大 ， 物 性 变化 的 影 
u y T 
图 5 变 物性 解 中 速度 分 量 、 等 温 线 之 间 的 比较 。 一 0 =0.0101， see 0,204035, ---6=08. 
Fig. 5 Comparison of velocity components and isotherms in variable property solutions. —— 60, =0.0101, =- 0, = 0.403 , 
--- 6 =08. 
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图 6 外 圆 环 上 传 热 系 数 随时 间 的 变化 。 


Fig. 6 The temporal evolution of heat transfer coefficient 


calculated on the outer wall. 


图 6 给 出 了 不 同 温差 条 件 下 kK,, 随时 间 的 变 


eqo 


45. 258 =0.0101 时 ， 物 性 变化 相对 比较 小 ， 此 时 的 


变 物性 解 与 常 物性 解 之 间 的 相对 误差 很 小 ，K。。 略 


大 于 常 物性 解 中 x,, 的 值 。 随 着 温差 比 的 增 大 ， 物 


eqo 


性 变化 的 影响 逐渐 增强 ， 当 6 =08 时 ， 变 物性 解 中 


Ka 比 相同 Ra 和 Pr 条 件 下 常 物性 解 中 的 &,, 大 


eqo 


15%。 这 个 定量 的 结果 与 空 腔 内 的 结果 比较 接近 。 
这 表明 : 一 方面 ， 变 物性 的 影响 对 设备 几何 结构 的 
变化 并 不 敏感 ， 另 一 方面 ， 在 数值 计算 中 ， 通 过 忽 


略 流体 热 物性 的 变化 而 得 到 的 常 物 怕 


E 解 低估 了 对 流 


换 热 系 统 的 传 热 性 能 ， 而 且 温 差 越 大 ， 由 热 物性 变 


化 引起 的 偏差 也 越 大 。 
3.3 不 同 Ra 条 件 下 变 物 性 的 影响 


eL 


—11— Const 
5 ®— 0=0.0101 
4— 070.403 
了 一 0=0.8 


eqo 


100000 
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图 7 不 同 瑞 利 数 条 件 下 外 圆 环 上 传 热 系 数 的 比较 。 
Fig. 7 The heat transfer coefficient obtained under different 


Rayleigh number conditions. 


为 了 进一步 分 析 变 物性 的 影响 ， 本 文采 用 变 物 性 
LBFS 模拟 了 不 同 Ra 条 件 下 同心 环 型 腔 内 由 温差 驱 
动 的 自然 对 流 ， 计 算 了 稳 态 阶段 外 圆 环 上 的 传 热 系 


数 K，。 图 7 给 出 了 不 同 Ra 下 变 物 性 解 中 的 K 


eqo eqo ? 


zi 
z 


由 图 可 见 ， 随 着 Ra WAK, ac 逐渐 增 大 ， 流 体 热 


物性 的 变化 对 x, 的 影响 逐渐 增强 。 换 言 之 ， 在 相 


eg0 


同 温差 条 件 下 ，Ra 越 大 ， 变 物性 对 传 热 的 影响 也 
越 大 。 


4 结 论 


本 文 提出 一 种 基于 流体 物性 变化 的 格子 - 玻 尔 效 
曼 通 量 求解 器 〈 变 物性 LBFS)， 该 求解 器 保留 了 传 
统 LBM 的 优势 ， 同 时 克服 了 传统 LBM 在 网 格 划 
分 ， 边 界 条 件 的 处 理 等 方面 的 缺陷 。 有 关 同 心 环 型 
腔 内 由 温差 驱动 的 自然 对 流 的 数值 研究 结果 表明 : 
通常 关注 的 常 物性 解 低估 了 换 热 设备 的 传 热 速率 ， 
尤其 是 在 温差 较 大 ， 瑞 利 数 也 比较 大 的 情况 下 。 
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